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A.JJLXJJUSJE 

With the advent of microprocessors) the DFT computation 
for spectral analysis has become very easy and economical ¥ But 
to carry out spectral analysis for real-time applications the 
choice of appropriate algorithm is imminent Recursive way of 
computing DFT 1 is the most suitable one for the said purpose* 
After briefly reviewing the various algorithms for computing 
DFT, this work discusses the implementation of a known recursive 
algor tihm on an Intel 8080 microprocessor based computer 
CECIL'S MICRONS) The advantages of this algorithm when compar-* 
ed with FFT algorithm are clearly brought out in this work by 
making a comparative study of both the methods Several examples 
are considered to demonstrate the utility of the algorithm in 
computing the DFT recursively A specific mention is made about 
the inability of the developed assmebly language program for 
carrying out spectral analysis of high frequency signals because 
of the short comings of the microprocessor used An economical 
and highly reliable real-time spectrum analyzer scheme is 
proposed for actual implementation In conclusion an outline 
for further work based on the concept of recursive algorithm 
is also suggested 
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CHAPTER 1 
INTRODUCTION 

In recent years digital signal processing techniques have 
gained great importance in many areas of science and technology 
This^mainly due to the fact of increasing availablifcy of medium 
and f^ge scale digital integrated circuits with their desirable 

A 

properties of small size, low power requirement, low cost, noise 
immunity and high reliability 

Spectral analysis 1 of discrete time signals is a very 
important area of application of the digital signal processing 
Spectral analysis, in general refers to the various alternative 
descriptions of a signal in terms of a linear combination of 
chosen basis set of representation Fourier series and Fourier 
integral forms of representation are of fundamental importance 
because the exponential functions are eigen functions of a 
Linear time in-variant (LTI) systems Spectral analysis has 
got many diversified applications The major areas in which 
spectral analysis is extensively used are Electronics and 
Communications Engineering, Mechanical Vibrations and Structural 
analysis and in general for time serios analysis in all 
endeavours of scientific and technological activity to detect, 
identify and characterize the natural and man made processes 
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which change in time or space 

For example a high fidelity amplifier should not distort 
signals in the range between 50Hz and 20K Hz Spectral analysis 
is useful for checking this Another example is the use of sound 
spectra to identify the person whose spoken voice has been 
recorded This is a speaker identification problem which is mainly 
useful for crime detection Transmission system's properties can 
also be described in the frequency domain so that the system can 
be designed for undistorted transmission of the frequencies 
present in the actual signal that is being transmitted Spectral 
analysis is used for detection of signals burried in noise 
In air traffic problems the frequency content of a plane's 
trajectory corresponds to the sharpness of bends in the plot of 
this trajectory on a radar screen The range of possible 
frequencies in this curve wall determine how close together in 
time we must take position samples for a digital signal processor 
which predicts the plane's future position A radio Astronomer 
may be able to deduce information in a planet's atmosphere on the 
basis of relative strengths of the planet's radiation in various 
frequency bands Mass spectrography and crystalo graphic analysis 
are among the many areas in which frequency analysis Is used 

Even for large amount ot data^ spectral analysis has become 
a considerably easier affair with the advent of digital signal 
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processing either by digital computer or dedicated hardware 
Further it has become very economic with the advent of ISI 
technology and microprocessors The use of such devices for 
the purpose of signal processing necessitate tho description 
of the analog signal from various transducers in the discrete and 
digital domain 

Thus the concept of Fourier representation of analog signals 
is extended to data sequences and digital signals The result 
is *Discrete Fourier Transform 1 (DFT) This form is nothing 
but discrete-time and discrete-frequency representation of 
Fourier transform 

With the advent of small, inexpensive mici^computers and 
with the developed algorithms for efficiently computing DFT, such 
as "fast Fourier Transform* (FFT), Winograd Fourier Transform 
(WFT), Fourier Transform computers using coordinate rotation 
digital computer (COHDIC) iterations and recursive ways of 
computing DFT, the utility and applications of DFT, has increased 
enourmously The applications include analysis of mechanical 
and structural vibrations, the extraction of Vital information 
from radar and sonar signals, analysis of statistical data, 
digital filtering and many others 

To carry out the spectral analysis in real time we have to 
compute the DFT as the data come and without any storage After 
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a survey of algorithms available to compute DFT in an efficient 
manner along with the constraint of ’real time spectral analysis’, 
’recursive way of computing DFT 1 is found tho most suitable 
So in this present work recursive algorithm is made use cf for 
finding out the DFT of a sot of sample points This type of 
spectral analysis is useful especially when the data are coning 
Sequentially and when there is sufficient tine gap between samples 
This algorithm has certain advantages as compared to computation 
of DFT using FFT algorithms A comparative study has been 
conducted and the results are presented m later Chapters 

The recursive algorithm for computing DFT has been Implemented 
on a micro-computer (ECIL’s MICRO-78) which is built around Intel 
8080 microprocessor This particular machine is selected mainly 
because of the availability of this microprocessor with a number 
of required interface, peripheral devices and support softwaro 

The DFT for different sizes of data, obtained from time 
sampling of a representative set of signals has been calculated 
on tho micro-computer with the holp of a recursive algorithm 
An implementation of the same algorithm on DEC— 10 computer 
yielded approximately the same results 

Some methods or algorithms available to compute DFT of 
a given set of data points and tho theoretical formulation of 
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the recursive algorithm will be dealt in Chapter 2 The 
details of microprocessor implementation of the recursive 
algorithm are presented in Chapter 3 Chapter 4 is devoted 
to illustrative examples of DFT calculation using the recursive 
algorithm Here as a particular application, DFT is utilized 
in the computation of power spectrum estimate which is a very 
useful characterization of stationary random signals 

Chapter 5 undertakes a comparative study of computational 
tirno and memory requirement aspects in the computation of DFT 
using a) recursive algorithm b) FtfT A brief report on the 
conclusions of the work and the scop of future work is mentioned 
in the last Chapter 
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CHAPTER 2 

THEORETICAL FORMULATIONS 

A number of algorithms are available to efficiently compute 
OFT of a set of data points We have to select one particular 
algorithm depending on many factors such as l) the sizo of data 
for which DFT is computed 2) the particular application in mind 

3) the spet-d of the processor at hand which computes the DFT# 

The following algorithms are available in literature for conpu» 
ting DFT 

1) DFT by use of FlfT algorithms 

2 ) The DFT algorithm of Winograd (WFT) 

3) DFT vising OORDIC (Coordinate dotation Digital Computer )> 
V?FT (very Fast Fourier Transform) 

4) Recursive algorithm for DFT 

5) State variable approach in calculating DFT (Recursive 
algorithm) 

In succeeding sections of this Chapter a brief mention of the 
various algorithms is made and finally the details of Recursive 
algorithm are given in Seotion (2 8)* 

2 1 THE FAST FOURIER TRANSFORM (F*T) 

The is simply an algorithm (a particular method of 
performing a series of computations) that can compute the DFT 
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nuch more rapidly than the direct evaluation of DPT 
DPT is defined as 

W"1 4)r 

s ^ 2k W tea) 

3*0,1, M-1 

vhero 

W * exp (-2TT1/M) 
i * A1 

The computation of DPT can be written in matrix form as follow* 
(for H « 4) 

Vf° W° W° W° 

VJ° W 1 W 2 W 3 . 

W° V/ 2 W 4 W 6 

W° w 3 w 6 W 9 

or more compactly 

X(n) * W* k Z(k) 

2 

If samples are complex numbors we need M complex multipli- 
cations andM(M-l) complex additions By suitably factorizing 
W matrix we reduce the multiplications The details may be found 
in Reference (l) Essentially for , the FFT algorithm 

is a procedure for factoring MXM matrix into y matrices each 
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(MXM) such that eadi of the factored, matrices has the special 
property of minimizing the number of complex multiplications 
and additions FBT requires My/ 2 complex multiplications 
and My complex addition, -where as the direct method requires 
M^ complex multiplications & M (M-l ) complex additions Thus 
large saving will be there in the computational time The 
above is one method of computing DFT efficiently 


2 2 FIT ALGORITHM DETAILS 

In the following section the fast Fourier transform (F^D^- 
al«orithn is briefly described This is considered worth 
describing because this brings out clearly some of the advantages 
in DFT computation by recursive way of computation 

Considering the DFT relationship 


X(n) * 


(N ® 2> Y ) 


Nrl 

k=0 


^Ck)w** 


(2 4) 


n 5:2 0 jlj 

where W * e"^ 2 tr /W 


N-l 


For H = 4^ y - 2 and k and n are represented 2 bit binary 
numbers 

k B 0, 1, 2, 3 or k = Ck x , k 0 ) 

a 00, 01, 10, 11 



9 


n * 0, 1, 2} 3 or n » (n^, n^) 

a 00, 01, 10, 11 

A compact method of writing k and n is k = 2kj_ + k Q , 

n a 2n^ + n^, where k 0 , k^, and can take on values of 0 and 1 

only 


Using the above representation for H*4 the DFT can be 
written as 


X(ni,n ) m 


1 

i 

k 0 =o 


Jo 


x 0 & 


a* 


lc 0 )W 


fenj+n^teKj+Kg) 


<2 S) 


It should be observed that single summation in DFT 
original definition, equation (2 4) is now replaced by V 0»2) 
summations in equation (2 5), in oraer to enumerate all the 
bits of the binary representation of k 

Factorization of w p 

w &ni ♦ V (2k x + k 0 )^ (2n x +■ i^) 2k x w (2n^ + 


(W 


, 2n o k i 


) (W 


Cw 


(2n x + r^) k 0) 


2n„ki (2ni + n ) If® 
(W 0 1 MW 0 > 
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This is because, 'in.ki 4 rLk*. 

’ w * =(W ) *■ 1 


/ - a \ 

r 4 j 

Then the equation (2 5) can bo written in the forn 


*l k l 


, ^l k l 
= (1) - 1 


XCn^i^) 


I 


v° 


1 

t X 

k i=o 


w (2ni+n ° )k <3 


(2 6 ) 


Equation (2 6) represents the foundation of the FtfT algorithm 
To demonstrate this point, lot us consider each of the sunnationp 
of equation (2 6) individually 

If the summation in the brackets is rowritten 


1 ^ruki 

^l^ko* 83 *o (k i,k 0 ) W 

X x (0, 0) 0) + ^(1,0) W° 

%<(>, l) =^(0,1) 4- Xq(1, 1) W 

X x a, 1) = 3^(0, 0) + Xq(1> 0) W 2 

a^d, i) -SqCo, i) + *od, i) w 2 


(2 7) 


If those equations are written in matrix forn as 



J 


V°> Q) 
*1<0, l) 


1 0 W° 0 

0 10 w° 

' 

— , 

Sf 

o o 

o o 

H O 

^ ~ l 

! *i a » °> 


1 0 VJ 2 0 

1 

0) I 

1) 


0 10 V/ 2 

* - 

k 

*0(1, I)J 


Thus, the inner sunnati on of equation (2 6) specifies the 
first of the factored matrices for the example 

Similarly, if the outer summation of equation (2 6) is 
written as 


Xgfaio, n x ) = 


\ 


k o =0 


%( n o» lt o) W 


fenj+r^) k c 


and can ho written in^ matrix form as 3 

O 


p*» 'Ml 

> 

£ 

o 

** 

o 

V * 


Xq (Qj i) 

: j 

55 

31 <1, o) ; 


*2 ( 1 . 1 ’ 
m* -M 



1 

1 

0 

0 


w 

w 2 

0 

0 


0 

0 

1 

1 


0 

0 

w 1 

v 3 


(2 9) 


. 

%< 0 * °> 
* x ( 0 , l) 

- 1 

%< 1 , o* 

J 

i 

\ 

*L&. 1 ) 

0 W* ' 


(2 10 ) 


Thus the outer summation results in tho second of the factored 
matrix 


Erom equations (2t6), (29) ono can observe that 
X(ni, Hq ) * « 2 (n o» n l> 
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If the equations (2 7), (2 0), <2 11) are combined th* foil- 
owing set of equations will be obtained 

/ s ^ 2n n k i 

*l (n o> V = 1 *o(ki,ko) w ° 1 

*S<"d) Kl> = l *]> 0 ,ko> W (SQl + n ° >k ° <2 12) 

0 

X(n 1? i^) = ig (r^, r^) 


This set of equations (2 12) represent the original 
Cooley - Tukey algorithm of F£T 


The above equations are represented by the signal flow-* 
graph and is shown in Fig 2 1 

XqOc) X x (k) Xg(k) 



X(0) 

X(2) 

XCl) 

X<3) 


Fig* 2*1 


13 


It is observed that input data is in natural order and 
output is m scrambled form One has to interchange nodes 
(01) and (10) in each array to get output in natural order 
and input data then becomes scrambled one In FPT algorithm 
there should be one sub-routine which calculates the bit 
reversal of a numbers (unscrambling) This difficulty is not 
there in recursive algorithm to be dealt in later sections 
After scrambling the input data the signal graph will look 
like as shown in (Fig 2 2) 


XoCk) X x (k) X^(k) 



Pig 2 2 


The two forms of the algorithms are referred by the term 
decimation in time Th 0 terminology arises from the fact that 
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alternate derivations of the algorithm are structured to 
appeal to the concept of sample rate reduction of throwing 
away samples and thus the term decimation in time 

Sande-Tukey developed similar^?* algorithm In contrast 
to Caoiey-T ttkey approach here components of n are seperated 
instead of the components of k Finally here the following 
equations are obtained 


(ho» 


k 0 ) = 


Jo 


*0 & I1 




w 


2n o k l 


V 


Vo 


■*2^o> ^1^ 


Xo 




V W 


2n l' k o 


XGi^no) = XgCn,,,!^) 

This form of the algorithm is known by the term deotuoation 

atfjdt 

in frequency jjbhe reasoning for the terminology is analogous to 
that for decimation in time Other details can be found In 
ftef (l ) • 
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2 3 FOURIER TRANSFORM COMPUTERS USING CORDIC ITERATIONS 

The basic Discrete Fourier Transform (DPT) is defined as 


M~1 

w, = y Zrr (W) 

3 K=0 * 


JK 


j - 0 ? 1, 


M— 1 


where 


V/ * exp (-32 it i/M) 

i » (-1)* 


(2 14) 


Zq> Zjl %m- 1 are samples of the waveform 

In this method of computing DFT the transform may be viewed 
as a process of rotating the vector a K (real and imaginary 
components) by* the angle (2 tt jK/M), followed by a sum over all 
K for each j The CORDIC vector Rotation method of Voider 
Ref (2), is used for the computation of DFT This appears to 
be attractive as a vector can be rotated in the time of a 
single multiply using this technique 

In addition, trigonometric coefficients as distinct from 
the CORDIC constants are not needed in a CORDIC rotation 
procedure so the expense of generating, storing and retrieving 
these coefficients is avoided 
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2 4 GORDIC TECHNIQUE 


Consider the vector X± = (x,y) to be rotated by angle 
= 2 it jK/M Let a new vector ^.i+i be obtained from Xj, by 

the process 


=H+1 = xi + y* 8" 1 

*1+1 = - *i s' 1 

Then if 

% = x i 2 + y ± 2 

= Tan^Cy^/xj^ 


(2 15) 


(2 16) 


the result will be that the magnitude and the angle of the 
new vector 


^i+1 “ kcL+1 > 

will be 

%+l = Hill+2 -21 1* 


6 i+l = 0 1 " tm' 1 (2' 1 ) 

Now to rotate X by an arbitrary angle 0 between (+ n/2 ) 
to an accuracy of 1 bit in n bits, the following procedure 
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way be used 

Initialize 1=0 

Zl - - * 

Iterate n times 
23 Sgn(Z^) 

^i+1 K \ + ^1^1 s”** 

71+1 55 ^1 * a i *1 2 "* i 

z i+l ~ %± ~ \ tan - ^ (2~* ) 

i » i + 1 


These iterations are composed of addition, subtraction,. shifting 
(multiplication by 2“*) and table look up (tan"*^ (2”*' ) ) 

2 5 APPLICATION TO DFT 

The application of CORDIC technique to DPT is as follows 

One rotation will be required for each combination of 
input and output for N 2 total revolutions associated with N 2 
multiplications in the computation of DFT The real and 
imaginary" components of the input vector are taken as the 
Initial values of *i and y^* The initial value of is 

then set to 2 irjk/M The CORDIC iteration is performed n 
times to produce n bits of precision and th resulting x and y 
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values are added to the sum that will become Wj after M repe- 
titions of this proooss The process is repeated for each fre- 
quency component W-j We also note that the resulting spectrum 
is scaled by a constant factor k where 

n— 1 ~ 2 ± i 

n (1 + 2 r is (8 19) 

i=0 

A very low-cost Fourier transform computer can be cons- 
tructed to perform DFT as described above It will be most 
useful for those oases of modest input data rate, where only 
a few high resolution fr quency samples are required for the 
output in this case considerable storage of trigonometric 
coefficients , unnecessary frequency coefficients and buffers 
can be saved as compared with the EFT approach 

The basic idea of application of CORDIC technique to DFT 
computation has recently been reviewed by ALVIN 1} BBSFAlN, 

Ref (3) who has proved that by suitable modification of CORDIC, 
DFT computation will involve no multiplication The following 
description gives an idea of DESPAIN*s very Fast Fourier 
Transform (VFFT) using CORDIC techniques 

2 0 VERY FAST FOURIER TRANSFORM ALGORITHM 

This algorithm provides procedure for computation of 
DFT that Involves no multiplications As DFT can be viewed 
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as a sum of rotated complex vectors) GORDIG techniques can be 
employed to accomplish the rotations It was also shown in 
the previous section (2 4) that an overall constant change of 
magnitude of the rotated vector docs not disturb the fundamental 
nature of the DFT but acts as a constant gain The acceptance 
of an arbitrary gain factor simplified tho CORDIC algorithms 
considerably In this particular algorithm (VFFT) it is visua- 
lized that all the angles can also be rotated by a fixed) 
constant angle Then the gam constant is complex and, if 
necessary, can be compensated for, in a final correction after 
the total DFT is accomplished As a result, multiplications 
can bo eliminated and the number of arithmetic operations 
required in the calculation of DFT can be considerably fewer 
than by other methods This results in considerable saving in 
hardware Because the algorithms require no multiplications 
and require only N storage locations (the so called ^nplace” 
property) they are of great advantage in microprocessor 
programs for signal processing tasks The details of the 
algorithm can bo found in Ref (3) 

In the following paragraphs very fast Fcurier transform 
(VFFT) algorithm is briefly described 

I 

DFT is defined as 

yk 

Bjr W N r « 0,1,2 N-l (2 20) i 

t 


N-a. 


r 

[ 
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where 

W N e exp (—2 j/N) 

5 “ /-I 

Here it is assumed that N has at least one factor equal 
to 16 and the DFT should he calculated with as little effort 
as possible The radix-16 algorithm can be easily derived 
from equation (2 20) Given N ~ 16T, let r » s + Tt and 
K = 1 + 16m 

Now eqn (2 20) becomes 
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Thus 


l 5 It SI Sm 

A S + tT * , L W 16 W 16T I B l+16m W T 


1=0 


m=0 


(2 21 ) 


Now let 


SI TjtI 

B*„ , = i B_ 


Sm 


'8,1 - "161 m t 0 ”l+16m W T 


(2 22 ) 


and (2 21) becomes 


A S+tT * X B« fn 


1*0 


Sl 16 


(2 23) 


This equation (2 23) is the DPT for N a 16 while (2 22) is 
the DPT for N *, T multiplied by the factor This procedure 

can be applied to eqn (2 22) until all possible 16-point 
DPT f S have been realized (i e T no longer has 16 as a factor) 

The equation (2 23) can be manipulated and carl be conver- 
ted to radix-4 representation And finally the 4-point DPT 
is convex ted to 2 point DPT after repeating the same mathe- 
matical manipulation that yielded 4-pomt DPT ffcoin 16-point 
DPT Fig 2 3 is a block diagram illustrating ’recursive radix" 
representation of a 16-point DFT In the above derivation* 
the only term which involves more than simple addition or 


? 


i 
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subtraction is the vector rotation operator It can be 

converted to a form in which only a minimal number of additons, 
subtractions , and shift (multiply by 2^) operations occur 

Determination of *16 * 

Table I illustrates the vector rotations of W^g as a 
function of i The method used to calculate is to 

employ rotations of the form 

x* «■ fxCx) * f a (y) 

y r «• fiCy) t t 2 <x) 

where A =* x + /* y 

and f^ and f 3 are restricted so as to involve only additions, 
subtractions and multiplication by 2* (shift t g ) 

The resulting rotation ©- is 

&{.«= + tan"' 1 - ( i. )/f^ ( )) & 24) 

while the magnitude of A is multiplied by G » 

G = /ff( ) + f| ( ) 

For the CORDIG case ()=()* = ^ ^ 2 ~ 1 and iterations 

for i a o, 1, 2, n-1 are employed to rotate by a desired 
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angle with n bits of precision in the result This requires 
n additons and n subtractions Another approach Is to imple- 
ment the four rotations 8i> % and shown in Table I 

and compensate later for the bias 

Rotation by + CM 

To obtain the rotation by + ft then form equation (2 24) 

£6 =? tan“^( P/ a) or fl/a « tan(fr) 

One has to choose and p to produce 0 to the desired degree of 
precision, and such that the corresponding equations 

x* «■ f-^x) + f 2 (y) 

y« «■ f^Cy) + f 2 &) 

contain only a minimum number of additions and subtractions 
Solutions for various levels of precision were generated by 
use of the following algorithm 

Algorithm 

°Find integers 
bit in N bits 1 
Begin a 

e 


a p such that a / p = R to within + £ 


a 


N-l 


8 **) 


while ( a 
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Begin a + a + 1 

$ «* Integer part (a R ~ e) 

while (f3<aR+e+ l) 

Begin 

Error ♦ | ( 0 / a ) - Rj 

If (Error ^ e ), Display(a , Error) 

P «• B +1 

End, 

APPLICATION 

A Radis - 16 Pipeline Cascade 

To illustrate the application of the above principle 
in an F FT processor, a 16 point DFT cascade circuit with a 
precision of 1 bit in 16 bits (dynamic range = 90 dB) will 

be derived This circuit could be employed within processors 
that provide DFT of any size which includes 16 as a factor 
The basic form is shown In Fig 2 3 

The 2-point DFT is of the form 

} rK 

^ % W 2 

Thus Aq =5 Bq + Bi 
and A]_ = Bq *» B]_ 
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A digital circuit to implement those operations is 
shown in Fig 2 4 

Thu rotation W 4 of Fig 2 1 can he expressed as 


uf 


VL 


(exp(- gir i/4)) 


uf 


ir j 

= exp ( Uf ) 

2 


Thus 


W 


uf 


or 


w; 


W, 


V 


w 


w 


4 


(~i) uf mod 4 
1 


-i 


-1 
+ i 
1 


A digital circuit to implement this equation is shown in 
Fig 2 5 
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The last operation of Fig 2 1 is Let i = uf $ then 

one can use the representation of Table I for W^g The rota- 
tions &2 are implemented with circuit of Fig 2 £ as 
explained above rotation by may be implemented by referon.ce 
to Table IX Ca) To obtain at least 16-bit precision, eight 
adders will be needed to realize th e ratio 

(f3 /ty - (138/309 ) decimal = 


(0 010000000 / 0 0100110101) binary 

then x f and y r are given by 

x 1 «• 0 01000110 x + 0 lOOHOlOy) 

y l «• o lOOHOlOly i o OlOOOOOOOx 

or 

x' = «x + X2" 2 ) - he* X 2 -2 ) 2~ S ) X S’ 8 + y 2“® 

yl = ((y + y2 -S ) - (y + y 2 -2 ) 2 -5 )-62“ 8 ± Xg* 2 

**2 

and if we define tx = x * x2 

ty = y + y2~ 2 
then the equations become 

X ( = (tx - tx 2*“ 5 - X 2~ 8 ) 2~ 1 + y2~ 2 

y 1 * (ty - ty2"* 5 - y2~ 8 ) 2 _1 £ x2~ 2 
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and it Is apparent that eight adders all needed Fig 2*£ 
illustrates the ckt rotation by is accomplished in a 
similar manner using Table I & II (b) 

Here (f^) „ ( VO) C(~2"* 0 )) 

x 1 (x + x 2 ~ q ) + y 2“ 2 
y 1 * (y + y s~ 2 ) + x 2~ 2 

x 11 ♦ x 1 + y 1 2*" 10 

y 11 * y 1 + x 1 2“ 10 

The resulting circuit is shown in Fig 2 7 complete 
pipeline circuit for the 1^16 cascade is shown in Fig 2 S' 

Thus the VFFT algorithm is more suitable to hardware implement- 
ion For hardware Implementation of VFFT algorithm, making use 
of recently evolved CORDIC arithmetic processor chip is novel 
idea* A monolithic processor developed by Gene 1 Haviland and 
A TUSZYNSKI, computers produo ts , quotients and several 
transcendental functions A bulk of CMOS technology with 
conservative layout rules is used for the sake of high relia- 
bility, low power consumption and good cycle speed The 
details of this chip may be found in Ref (4) 

The details of Winograd algorithm to compute DFT are 


availably in Ref (5) * 
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2 RECURSIVE ALGORITHM 

There are many ways of computing DPT recursively The 
particular method will he selected depending on the applica- 
tion The details of state variable approach to recursively 
compute DFT can be found in Ref (6) 

In some applications it is convenient to segment the 
complete signal and do the spectral analysis segment wise 
This is especially so in case of a sophisticated signal f (t) 
of long duration, for example, a voice record In principle 
we can uniquely determine f (t) in terms of its transform 
“p (W), however any numerical processing directly involving 
F(W) is unrealistic because F(W) is too complex to be recovered 
with sufficient accuracy from its samples To avoid this 
we propose to describe the spectral properties of f(t) 
inform. of the running Z-transform or recursive DFT formula 

The DFS (Discrete Fourier Series) of N numbers can be 
obtained from the running Z Transform The following formu- 
lation will yield a method to compute DFS in a recursive way 

Recursive DFS 

Given a sequence f (n) and an integer N, we form the 
running Z transform 

N-l . 

0 (n,Z) = f(n-K) Z * 
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of f(n) For a fixed 'h* 0 (n,Z) is the Z transform m the 
variable K of the segment f (n - k) of f (n) shown in. Fig 2 10 



Fig 2*10 

The function 0(n, Z) satisfies the first order recursion 
equation 

0 (n,Z) - IF 1 0<n - 1, a) (2 26) 

= fCn) - Z" N f Cn-N) 

This can be proved as follows 

N-l v. 

0 <n, z) ss f (n*k) z 
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N-l 


0(n-l,Z) = l f |n-(k+l)^** K substituting k+1 = r 

k=0 “ . . 

k= r ~ 1 


0 (n~l , a ) « | 


^ ~(r-l) 

f(n-r) Z 


N 


» Z 


0(r-l,-E) 


f(n-r) £ 


-r 


N -r 

l f (n-r) 3L 


r=l 


+t 

0(n,Z) - Z 0 (n-l, Z) 

N-l N 

= t l £ (n-k) Z I - t £ f(n-r) Z r J 

k" - 1 3?"“1 


= f(n) - f (n-N) z“ N 

Therefore, the resulting recursion becomes 
0 (n,Z) -*" 1 0(n-l,Z) 

= f (n) - Z~ N f (n-N) 

N-l 

0 (n.Z) = l f(n-k) Z K 
k=0 


if yre substitute Z = ) - w 
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0 (n, W"“) =: l f (n-k) W (2 27) 

k=0 

Equation (2 27) shows that for a fixed n, the DFS of the N 
numbers f(n*4c) equals 0 (n, V/" 01 ) 

The running DFS 0 (njW"®) of f (n) satisfies the recursion 
equation 0(n,W“ tn )^0(n-l, vf®) = f(n)-f(n-N) (2 27) 

This follows from the previous recursion formulation 
(eqn 2 26) with Z = V/’ fIn and using the fact V/ N = 1 

From this recursive equation DFS analyzer is construe- 
t ed as follows 

The recursion equation (2 27) defines a discrete system 
with input f(n), output 0(n, W“ ra ), and system function 


S(Z,m) = — ; a (2 28) 

1 - vP 1 z 1 

The system is shown in Fig 2*11 



t 


Fig 2 11 




The system consists of a shift register with output 
f(n-N)* one delay element and multiplier The shift register 
can be omitted if we have direct access not only to f (n) hut 
to f (n-N) Connecting N such systems in parallel we obtain 
the real-time spectrum analyzer shown in Fig 2 12 



F m =3 DFS Terminals 


Fig* 2*12 
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If the Input to the analyzer i s a sequence f (n), then 

the outputs at terminals F m equal the DFS coefficients tfCnjW** 111 ) 

of the N numbers f(n) f(n-N+l) 


Example 

Let us say the sequence for which DFS is required to 
be computed is 


f(-5), f(-4), f(-£>, fC-l), f(0), 

f(l), f(2), f(3), f(4), f(5) 


Let us say we want to do the spectral analysis for each 
segment having 4 points Binning Z-transform enables us to 
segment it as follows 

segment (3) 
up 3 


l 



(-5),f(~4),f (~3)»f (-2),jr<-l),f (0),f (l),f C2),f (3),f(4) 



Segment (4) 


If the DFS of segment (4) is to be computed by making 
use of recurive formula it can be found from DFS of s egment (3 ) 
and Ofte^new and old samples t 

t 

f 


\ 
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is a numerical example consider the sequence 

2 ^ 5 j j 3, 7 y 1 ^ 2^ 1, -5 

f(-4),f(-3), f(-2), f (-1), f(0), f(l), f( 2), f(3), f (4) 

N = Segment length = 4 

To compute, let us say, 0 (4, W ) for m - 0, 1, we 
can illustrate the u e of recursion as follows 

N-l w. 

0 Cn, VT 111 ) = l f (n - kJ W 
k=0 

32 */N 

N = 4 , N-l = 3 V - e 

V- . 0 J 2 v / 4 = e 3^/ 2 = j 

3 

0 (4,w”°) » j> f (4 - k) W° 
k=0 

» f (4) + f(3) + f(2) + f (1) 

£=-5 +1 + 2 + 1 - -1 

This is actually the segment using samples f(4), f(3), 
f(2), f(l) 

»*4 0 ( 4 , w*°; by definition has got the value -1 

Similarly by definition the value of 0 (4, W**" L ) is *1 *■ 
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These two values can be computed recursively from 0(3, W"°) 
and 0 (3, W If we calculate these two values we get 

them as 11 and respectively (can be calculated from 
definition of 0(2,Vf^) and 0(2, V/** 1 ) 

Recursive formula states that 
0 (n^) - W m 0 (n-1, W“ m ) - f (n) - f(n-N) 

0 (4,W ) can be obtained from 0(3,W ) by putting n = 4, 

m=0 m the above equation 

0(4,vf°) s W° 0(3,vT°) + f (4) - f( 0 ) 

- 11 + (-5) - 7 
= - 1 

Similarly 0(4, W” 1 ) = W 1 0 fejN" 1 ) + f (4) - f(0) 

_ j («*j5 ) + (-5) — 7 
= »7 

Observations » Here one has to fix the segment length The 

sequence length for which DPS is computed segment wise, needs 

be 

to /fixed, if the record of the data before a reference time 
is not available In practice a sequence of zeros of length 
M~l is augmented before the data sequence so that only the first 
sample is needed to start with 

An equivalent recursive formula can be derived directly 
from the definition of DPT This should be so because of 
the fact that from Z-Transform one can get DPT with proper 
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transformation The recursive formula derived directly 
from the definition of DFT is implemented on Micro computer 
in this work The advantages of recursive computation of 
DFT are also enumerated in the next section along with its 
formulation This recursive Way of computing DFT is useful 
especially in situations wherein data for which the DFT is 
computed are received sequentially and computation of the DFT 
can not he completed until the final datum is received This 
is especially suitable n there is sufficient time gap 
(computational time) between samples By computing the DFT 
recursively, only a small number of operations need be made 
after the receipt of the final datum 

In the following section the recursive formula is 
derived# This particular form is used in doing the micropro- 
cessor based spectral analysis m this work In the next 
section advantages of this algorithm when compared with F3fT 
are enumerated 

2 8 Derivation of Recursive Formula* As we know for a sequence 
z 0* z l* of sampled waveform data (either real or 

complex ), the DFT is defined as 

M-l 

Wj = y Z, exp (-2 rr ijk/M) (2 29) 

J k=0 ^ 


i s/-l 


42 


for j = 0,1, M-l An assumption implicit in the defi- 
nition of W-j, when applied to sampled waveform data, is 
that the samples Z 0 , Zi Z^, are obtained at equally 

spaced time intervals Tq 


The Moving -window Discrete Fourier Transform (MWDFT) 


Let Zq, Zi , Z m be an infinite sequence of 

sampled waveform data, where Z]£ is the (possibly complex) 
value of the signal Z(t) obtained at time t = K T 0 Suppose 
that for n = 0, the DFT of Zji, Zj^i, Zn+M-1 i s 

to be computed, this sequence of DFT*s is referred to as the 
moving -window DFT (MWDFT) 

For n = 0, 1, MWDFT can be defined as W^ n as 




k=0 


for J » 0,1, 


Z k+n exp (-2 * ijk/M) 

i * /-i 

M-l 


M-l 


Wj, n+1 = J 0 Wl ■**<-«*«»» 

Substituting K+l = n 

K s* K 1 -! 


* 
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w 3,oH 

= exp{* 


M 


= i4i 
UAi) 

M 


Z K i +n exp — (K 1 -!) 

K il 2 h\n exp ( — sr — > 


3nlj 

e M t Jf Zjl exp( ~~ )] + Z(n) - Z(n) 1 

K — 1 * + n « 


SjJiA M-l -2 ijk 1 

3 e kio Zk1+o exp ( ~ ir_) ' z<n) + w 

Ther afore the desired recursive formula is 
^j,n+l *^ w Jn + z M+n " Z(n) 1 exp ( 

Wjo = VJj 

for j *= 0,1, M-l, where is as defined in equation <2 29) 

It can be seen that for each n, the calculation of 
~ Z n ) is independent of the index j of the Fourier 
coefficients Thus, the total number of operations (an opera- 
tion is defined as the performance of a complex additon followed by- 
complex multiplication) for each n is just M, neglecting the 
one-time calculation of Z]^ n - Z- n The required computational 
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speed is therefore such that the M operations must be 
accomplished in T 0 seconds for a real-time application Thus 
we have to choose highest frequency of the signal appropriately 
keeping in view of the computation speed of the particular 
processor ^up Intel 8080 in the present work) at hand In some 
situations hot all of the Fourier coefficients are of interest, 
in such oases, the total number of operations to be performed 
is decreased and time per operation can be increased accordingly 

If one looks at the way in which DFT is computed by making 
use of recursive algorithm the following advantages ~an be 
attributed to this type of computation 

Advantages We can process the data as they come We need 
not wait until the final sample comes By computing the DFT 
recursively* only a small number of operations need be done 
after the receipt of the final datum Even if the recursive 
method of computing the DFT requires many more operations than 
would the FFT, the required time measured from when Zq is 
obtained, might be insignificantly different (That Is to say 
that the number of operations needed to compute tfca DFT by 
s4pe method is not an accurate measure of that method*s com- 
putational efficiency) Here, in this method we are saving 
time as there is no question of waiting for the arrival of last 
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sample and we are pro cessing the data as they came in sequence 
If ail the samples are stored as in the case of FTT the time 
required to retrieve the samples from memory area should also 
be included in the total time computation to have a realistic 
idea But in recursive way of computation there is no question 
of storage of samples and accordingly there is no such diffi- 
culty of frequent memory references Even the memory area 
required is small when compared with that required in PET 
All these statements are supported by quantitative results in 
Chapter 5 

An added advantage of this method is that the total 
number of points need not be a power of 2 

2 9 COMPUTATIONAL PROCEDURE - AH OUTU 

MWDFT is tsie’ fr 

V = io ^ n axp 

J a 0, 1, M - 1 

MWDFT for n equals zero is nothing but ordinary DFT 

Wj 0 = Wj 

The recursive formula states that 

, 2 *i;K 

w i»n+l = ' w jn + z M+n " z n J «M < ~J~ } 


i = / -1 
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Wjo * can be computed from the equation 

M»*l 

Wj - l Zfc exp ( ~ 2 .. y . J . k .) 1-0,1, M-l 

by simply starting with M-l zeros proceding Z 0 , the first 
sample in the data sequence The computational procedure is 
described below for M * 4 ( 4 point BFT) 


«* « 

Let the data sequence be Z^, Z^, Z^ y^start with 3(=M-l) 
z«yos preceding Z^ So the new sequence may be designated as 


J 3> 


Z4. 


2 55 z 6 


(Zq j Z X Zq ) 


z 3 > 24) * Z2M.2 correspond to the actual data sequece 

z o> 

First we compute W-jq - W-j for j = 0,1, M-l (=3) by 
making use of the formula 

3 -i2n 

w 30 = w d = Jo *k e ^~ ik 

It can be observed that while computing for 3=0, « 3 

(saM-l) only Z^ ml CZ3) has got non zero value* So the simpli- 
fication in computation at this step result This will give 



Then taking the recursive formula 
n+1 = * W 3n + z M+n * z n * 

<x>ui 

y first n=0 is substituted and get Wj^ for 3=0, M-l 
by making use of W-jq for 3=0, M-l and the next sample 
in the sequence i e , ZM+n ( = Z4+0) Finally when rt=2 (rr=M-2 ) 
we get final set of coefficients from W^3 which makes use 
of Wj2 and Z 2M-2 l ast datum The final set of coefficients 
which we compute after the arrival of last datum are exactly 
equal to the coefficients that are obtained by making use of 
the definition Recursive algorithm is implemented on 
MICRO-78 computer and results are obtained 


27 j 
M 
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CHAPTER 3 

microprocessor implementation details 


A recent and fascinating innovation in the digital 
electronics field is the development of the microprocessor 
using large-scale integrated technology The microprocessor^ 
combined with the memory and input/output devices, forms 
a micro -computer whose cost is competitive with conventional 
random logic in an ever-increasing number of applications 
The micro-computer is not only finding use in snail systems, 
where it nay replace a mini -computer, but it is also opening 
vast new areas of applications where larger machines are not 
economically feasible 

The microprocessor can be used for real-time applications 
and dedicated systems can be d veloped making use of single 
board computer (SBC) and related transducers With SBC and 
transducers, A/D and D/A converters and support-hardware the 
spectrum analysis can bo done in real time This type of 
spectrum analyzer will be very cheap to realize using micro- 
processors, especially with the advent of fast and cheap 
processors which are easily available in the markot This 
realization will have the advantages of low cost and high 
reliability when compared with spectrum analyzers realized with 
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combinational and sequential logic circuitry 

In this Chapter the details of Implementation of the 
recursive algorithm on micro-computer are given The assembly 
language listing of th program can bo found m Appendix A 

Performance limitations (spe d,shorter word length, 
United addressing modes, fewer internal registers) limit the 
competitiveness of micro -computers So while soleeting a 

microprocessor for a particular application we have to make 
suitability study of the processor at hand, and then only one 
particular type of processor should be selocted from the few 
which are readily available 

3 1 CHOICE OP THE MICROPROCESSOR 

Motorola 6800, Intol 8080 and Fairchild F-8 are tho few 
readily available processors For this project work Intel 8080 
processor is selocted The choice of Intel 8080 was based on 
tho standard support software available with it An Intel 8080 
based ECU'S micro -computer MICRO-78 with the facility of a 
floppy disk, resident assembler and also with an efficient 
debugging aid (on-line Debugging system) is chosen for the 
DPT algorithm implementation Description of the MICRO-78 
and Intel 8080 processor for our application is given in 

Hef XUO. 
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In this computer the data -word length is 8 bits But 

to represent numbers with reasonable degree of accuracy it is 

required to have 16 bits This is so because 8 bits are to 
for 

be spared/fractional part and another 8 bits for integer part 
of the number So two registers (Register pair) are used while 
doing arithmetic operations 

3 2 IMPLEMENTATION DETAILS 

The DPT computation requires values of sinusoidal and 
Cosinusoidal function values These can be generated by making 
use of trigonometric function routines provided in Decimal 
Arithmetic Packages (DAP) in MICRO 78 But this is nob preferred 
because the computation of trigonometric values by making use 
of the above said routines will reduce the computational speed 
of the processor drastioally* It also increases the storage 
requirement Instead?all these values, depending on the number 
of points for which the DPT is Computed are stored in the form 
of Look up Tables in certain blocks of MICRO-78 memory area 
These values are made use of while computing DPT The placing 
of these valued in the overall memory allocation is'&own in the 
block diagram of memory structure for the recursive algorithm 
in Fig « 3*1 



012000 

012075 

014475 

100 

/too 


Subroutines 
MULT, TWO >T WOT 


Mam program 

ABt Wj 0 for j»0 1, M 1 

AB III Wj n , j & 0,1 M 1 

and for each 
n*0 ,1 M 1 


Data points and 
look up tables 
(size depends on M) 


OFT Coefficients 
(size depends on M) 


FIG 31 RECURSIVE DFT PROGRAM MEMORY STRUCTURE 
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3 3 ROUTINES USED 

MULT This Is a 16 x 16 multiplication routine 
While writing this program every care has been taken to see 
that the total time for 16 x 16 multiplication is the least 
Multi pli cant is 8 bit integer and 8 bit fractional part 
Similarly 16 bits are allocated for multiplier The result is 
32 bits 2 bytes integer and 2 bytes fraction 

But it is truncated to 1 byte integer and 1 byte fraction 
for simplicity as it is not going to affect the accuracy much 
This routine is called in the main program as and -when the 
need arises 

TWO* This is a routine which will compute the two ! s complement 
of a 16 bit number This sub-routine is made use of in the 
subtraction of numbers 

is 

TWOT t This/also a routine which is useful in the subtraction 
of two 16 bit numbers This will introduce the sign of the 
resultant number and stores it somewhere for further use 

3 4 PROCEDURE FOR COMPUTATION OF FREQUENCY COMPONENTS 

Here the procedure which is followed to compute Fourier 
Components is outlined In general the frequency components 
may be complex The use of complex exponential routine already 

UH 

ictNo A 
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available in DAP routines of MICRO-78 is eliminated for the 
same reasons as outlined earlier (computational speed aspect) 

The complex Fourier components are spilt into real and imaginary 
and they are computed seperatcly This will reduce the overall 
computational time and enable th program to work for real time 
applications 

The DFT is defined as 


Mrl 

W< a l 

J k-0 


% exp ( 


~2ir ijk 
M 


i * /-l 


J * 0,1, M-l 

Wj n the MWJQFT is defined as 


WJ 


n 



k=0 


z k+n ex P ^iok/M) 


J as 0,1, M-l 

03ie final recursive formula says that 

2 n ij 

Wj, n + 1 * t W^ n + Z M+n -Z n )e M 
4 = 0,1, M-l 1 = ✓ 


(3*1) 


(3 2) 


(3*3) 


jjS . tep , X 

Computation of W^ n for n = 0 and j * 0,1, M-l, 
we get a sequence Wqq, %)* %) %-l> 0 in ' ttx3,s S ^ 0P 

i 
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Each number is a complex number Ihey are computed as follows; 

W j0 = Wj 

M-l 

AS * E Z K (~2*tok/M) 

k=0 

t s 0, 1, M-l 

and because of the fact Zq» z^ 2 becomes zeros when we 

start with M-l zeros in the beginning of the data sequence, 


R( V = W 3 REAL PARI = XtCOS x(M-D) (3 4) 

ijn (Wj) = Imaginary Part of W< 

1 

= Zy.-L X Sill ( —i. xJxCM-1) ] (3 5) 

Wj = R(¥j) + j VWp (3 6) 


The assembly language program starting with lable and 
memory location 012164 will compute Wjo for j = 0,1, M-l 
(both Real and imaginary components) W^ 2 VIj,M-2, 

W-j, M-l are computed from the recursive formula for 


3-0,1, M-l 

w j, M ^ w dn + %*-n " z n* 



6 
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REALCWj, nfl) = [Real Z M+n - Z n ] OOS (%*. -j ) 

- (Imaginary »jn)si» (fi) (3 7 ) 

BtAO(Wj » *1) = IMAO(Wj n )OOS ~ i * lREAl(W }n ) + Z M+n - Z n 1 X 



The samples are assumed to be real m the implementation of the 
algorithm on MICRO-78 

REAL (W-j, n+*l) and IMAG (Wj, n+1) are computed accordingly 
for n=0 This gives values for 3 = 0 , 1 , M-l Simi- 

larly making uso of equations (3 7) and (3 8) Wj2» Wj 3 W^M-1 
are computed by substituting n=l,2, M-2 The final DFT 

coefficients are obtained when n=M-2 is substituted (Wj, M-l 
values for 0,1, M-l) 

The assembly program listing starting from memory location 
012471 and label ABlll gives these oomponets. 

Thus the final set of components W 0 ,m-1 1 W 1>M«1* 

W 2»M-1 > %-l»M-l w111 glve the DFT of samples Z 0» Z 1 < 

Z M-1 
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CHAPTER 4 
EXAMPLES 


This Chapter undertakes various examples, which are 
considered in this work to illustrate the DFT computation of 
a given set of data points making use of recursive algorithm 
The succeeding section, discusses briefly the selection consi* 
derations of various parameters in the DPT computation 

4 1 SELECTION OF D#C PARAMETERS 

It is assumed that the discrete-time signal is derived 
frow a continous time function which is a representation of 
a signal The major DPT parameters are as follows 


T 

F 


% 

N 


Increment between time samples (in seconds) 
sampling rate (in Hertz) = 1/T 

increment between frequency components (in Hertz )=* 
frequency resolution 

record length (in s econds > effective period of the 
time signal « 1/E 

folding frequency = fg/2 (in hertz) 

highest possible frequency in spectrum (in hertz) 

number of samples in the record 


In order to avoid aliasing, it is neqessary that 
f« > 2 fh 


(4 1) 
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This implies that T must be selected according to the 
relation 

* , l/2f h (4 2) 

The increment between spectral components (F), can be 
thought of as frequency resolution For a desired frequency 
resolution the minimum record lengtn tp must be selected 
according to 

t p ^ 1/F <4 3) 

The study of equations (4 2) and (4 3) leads to the 
conclusion that there is a trade off between the bandwidth of the 
system and the frequency resolution between successive samples 
To be able to handle signals of larger bandwidth, it is 
necessary that T be reduced For a given N, this would 
shorten the record length and, thus, decrease the frequency 
resolution Conversely, to increase the resolution, it is 

necessary to increase tp For a given N this would increase 

of 

T, which would restrict the bandwidth/operation 

The only way in which either the frequency resolution 
or the bandwidth of operation can be increased while holding 
the other constant is to increase the number of points N in 
a record length If fh ^ are both specified, then, N 
must satisfy 


N / Sf h /F 


(4 4 ) 
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In all tho examples considered in this work the number of ■points 
for which DFT is computed is selected based on the above dis- 
cussion 

4 2 VARIOUS EXAMPLES FOR DFT COMPUTATION 
Example 1 

100-point DFT analysis to a pure sinusoid of frequency 
500 Hz and maximum amplitude of 5 and 10 Hz frequency resolution 
The DFT spectrum (plot of jx(m)|vs m) is shown in Fig 4 1 
A peak, as desired, at the frequency of the sinusoid is observed* 

Example 2 

Spectral analysis of a rectangular window shown in 
Fig 4 2 

Here 128 points DFT is computed and the sp ctrum is 
plotted It is clear from the figure that the amplitude of 
the frequency components decrease as frequency Increases and 
the envelope has the usual theoretical form of SINCx)/* The 
plot of lx(m) I VB m Is shown In Hg 4 3 

Example 3 

Power spectrum estimate of a rectangular pulse shown in 
Fig 4.4 is obtained after computing the 9-point DFT by recursive 
algorithm The envelope of the power spectrum estimate is 




n ~ 


N 128 points 

Sampling interval 
0 5 second 


FIG 42 RECTANGULAR WINDOW FOR SPECTRAL 
ANALYSIS 



x(n) 



8 


n 


N 9 points 
Sampling interval 
1 second 


FIG 4 4 RECTANGULAR WINDOW FOR POWER 
spectrum ESTIMATION 
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shown in Fig 4 5 The penodogram estimate of the PSD is 
computed from following the formula, 

S^w) = i I x„ (w)| 

N 

where X N (W) is the DFT of N Points 
Example 4 

, i 

By passing a jpseduo V^te-Guassian noise sequence of 
uniform power spectral density through a low pass digital 
filter which is obtained from tho analog transfer function 
of the Butterworth filter* the magnitade square characteristic 
of the low pass filter is obtained The Butterworth filter 
characteristic chosen is 

G(P) = i 

1 + 1 4 142136P + P 2 

The 3dB cut off frequency of the filter is 50Hz The sampling 
rate of the system is 500 Hz The magnitude square character!— 
stic can be found in Fig 4 6 The variations in the power spectr- 
um of pseudo— random noise sequence generated on the computer are 
reflected as variations in the filter characteristic obtained 





Xf m)j 



FIG 4 3 OFT SPECTRUM OF A RECTANGULAR 
WINDOW (HALF CYCLE) 






* * 
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CHAPTER 5 
COMPARATIVE STUDY 

The comparative study of the two algorithms to 
compute DPT of a set of given data points namely 1) PFT 
2) Recursive algorithm can be done at two fronts 
l) Computational time required to perform DFT 2) storage 
requirements 

6 1 COMPARISON OF COMPUTATIONAL TIME 

While comparing computational time required by the 
two algorithms it should first observed that FFT algorithm 
gives DFT for a set of given sample points where as 
recursive way of computing DFT takes sample points one by 
one and displays the DFT coefficients computed for data 
upto that instant To be specific, FFT algorithm has to 
perform I^logg M complex multiplications to get the M -DFT 

p 

coefficients where as M operations are required in DFT 
by recursive algorithm But to be realistic, one has to 
start counting time from the arrival of the first sample 
uptil the arrival of last sample for computational time 
calculation This is nothing but the signal acquisition time 
This depends on the particular sampling rate usqd 

In this work a logio state analyzer (Type HP 1610) 
has been connected to MICRO-78 computer and the computational 
time required to compute all the coefficients at a particular 
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instant of time is obtained The timings for different sizes 
of data are listed in the following table 


Table 5 1 


No 

of Points (Complex) 

Computational time 


100 

677 m 

sec 


128 

722 m 

sec 


256 

1415 m 

sec 


To get the final DFT coefficients these times have to 
multiplied by the data size (acually less than this is observod)^ 
because in the program one particular loop has to be computed 
recursively M times For M=sl28 the computational time is found 
to be 7225 m Sec , or 7 Sec , 225 m sec If we take the 
computational time required to perform the same operation using 
FFJE^t is found from Ref to be equal to u sec 500 m sec 

This is oust computational time If one adds the acquisition 
time which has to be added for comparison purposes $ the time 
requirod for computing DFT by recursive algorithm may be in- 
significantly different from that of FFT The term insigni- 
ficantly different 1 is used because one can compare the two 
methods only for certain sampling rates of the signal Especially 
if the sampling interval is more and the data size for whioh 
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algorithm also 1 in-place* notion is employed the saving in 
storage area is obtained When we compare FFT in-place 
algorithm to recursive algorithm there is no improvement in 
memory size But in thi PFT algorithm computational speed 
decreases So one has to select particular type of algorithm 
depending on practical constraint 
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CHAPTER 6 

CONCLUSIONS AND FUTURE WORK 
Concluding Remarks 

Assembly language program for Intel 8080 microprocessor 
based computer (SOIL'S MICRO-78) for the purpose of implementing 
recursive algorithm for the computation of DFT has been 
developed in present work This is useful in carrying out 
spectral analysis for real-time applications Same algorithm 
is also implemented on DEC-10 computer and results are 
obtained 

This type of computation is suitable when the data are 
coming sequentially and the sampling interval is sufficient 
enough to perform M operations for an M-point DFT Here, 
as DFT is computed from the definition there is no necessity 
for unscrambling of the output or scrambling of input data as 
required in FFT algorithm So there is no requirement of bit 
reversal sub-routines in the computer programe Even if the 
recursive method of computing the DFT requires mor G operations 
than would tho FFT, the required time measured from the time 
when the first sample arrives might bo insignificantly different, 
because of no waiting time and fewer memory references 
Another advantage here is that M, the number of points, 
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nood not bo a power of 2 

As the DPT algorithm is implemented in this work on 
an Intel 8080 microprocessor based micro -computer (EGIL l S 
MICRO-78) and because of tho computational time to perform 
M operations is of the order of few hundreds of m sec the 
spectral analysis in real time is restricted to input signals 
of bandwidth loss than few Hz This/^due to low clock frequency 
and 8 bit word size of the processor For representing numbers 
with reasonable accuracy atleast 16 bits are required, and so 
doublo precision has been employed This reduces the compu- 
tational speed to some extent The 16 x 16 software nulti- 
ptli cation routine which is developed and used in this work 
has taken approximately 1 26 m sec This can be reduced to 
approximately 350 micro-seconds if one dovolops hardware 
multiplier using TRW 1010 J (for 16 x 16 bit multiplication) 
chips » As the number of multiplications mainly determines the 
total computational time, with the processor used m this work, 
only low frequency signals may be handled with real time 
cos traint 

Some recently evolved algorithms for computing DFT 
efficiently are also presented This/with a view of tho 
possibility of efficient and economic realizations of spectrum 
analyzers in future f 
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Future Work i 

Thoro is scope for further work bas^d on the recursive 
way of computing DFT This is nainly so because of the 
increasing number of real-tim applications One can actually 
build such a real-time spectrum analyzer by burning the 
devoloped softwaro program, for th recursive implementation 
of the DFT, into an EPROM (erasable progranablo read-only- 
memory, chips like 1702A) and using it in the realisation 
In this, software multiplication routine can be replaced by a 
hardwaro multiplier built using TRW 1010 J chip 

Using this typo of 16 x 16 bit hardware multiplier 
and single board computer (SBC) kits having fast processors 
like Intel 8086 or Zilog 8000, appropriate interfacing hardware 
and a display unit, one can build an economic and highly reli- 
able real-timo spectrum analyzer 

An attempt can be made to extend the recursion concept 
of for DFT computation to compute other orthogonal transforms 
like Walsh -Ha damard transform (WHT) There may be some ad- 
vantages in terms of storage requirements in a computer if 
WHT is implemented recursively As the basis set in WHT 
involves only -1 and +1 there may not be an increase in 
computational time because of the increased number of opera- 
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tions in the recursive way of computing WHT In this t/P© 
of computations one should strike a compromise between 
memory size and computational speed The economy in terms 
of memory storage and recursion in memory references in computing 
VIHT by recursion may even override the loss in computational 
speed This is an exercise worth trying 
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